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~_" We analyze the recently released CoGeNT data with a focus on their time-dependent 

^^ properties. Using a variety of statistical tests, we confirm the presence of modulation in 

5^ the data, and find a significant component at high {E^^ > 1.5 keVee) energies. We find 

*^ that standard elastic WIMPs in a Maxwcllian halo do not provide a good description of 

^^ the modulation. We consider the possibility of non-standard halos, using halo independent 

techniques, and find a good agreement with the DAMA modulation for QNa ~ 0.3, but 

i-Q disfavoring interpretations with Qp^a = 0.5. The same techniques indicate that CDMS-Ge 

' should see an 0(1) modulation, and XENONIOO should have seen 10-30 events (based upon 

^ the modulation in the 1.5-3.1 kcVee range), unless Lcff is smaller than recent measurements. 

'~~' Models such as inelastic dark matter provide a good fit to the modulation, but not the 

CO spectrum. We note that tensions with XENON could be alleviated in such models if the 

f — peak is dominantly in April, when XENON data are not available due to noise. 
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I. INTRODUCTION 

The CoGeNT Collaboration has recently published results from the first fifteen months of data 
taking [U [2]. Since their first data release more than a year ago, they continue to observe an 
unexplained excess in the spectrum of nuclear recoil scattering rate and now claim an annual 
modulation of 2.8(7 [2j. This preliminary evidence for a modulation is an important step towards 
determining the nature of CoGeNT's unexplained spectrum and has been claimed to be evidence 
for a ~ 7 GeV dark matter [3j. In this work, we present a comprehensive statistical analysis of 
the CoGeNT data and show that the modulation spectrum is hard to achieve with a conventional 
light elastic WIMP in a standard Maxwellian halo. 

Direct detection experiments such as CoGeNT search for the scattering of dark matter off 
nuclei in ground-based detectors. The spectrum of nuclear recoil energies depends on the mass and 
scattering cross section of the dark matter, as well as its velocity distribution in the Galaxy. One 
of the most distinctive features of a dark matter signal is that it should modulate annually due 
to the motion of the Earth's rotation about the Sun [4| . In particular, the flux of dark matter as 
observed in the lab frame is larger in the summer, when the Earth is moving in the same direction 
as the Sun, than in the winter, when the Earth's motion is against that of the Sun [5j. For a 
Maxwell-Boltzmann velocity distribution, the flux peaks 152 days into the year. 

Observing an annual modulation in a potential signal is a crucial step in confirming its origin as 
dark matter. Direct detection experiments face the challenge of distinguishing dark matter nuclear 
recoils from a list of potential backgrounds. In most cases, experiments utilize a combination of 
ionization, scintillation, or phonon signals to separate out nuclear recoils [6]. But the possibility of 
contamination in the nuclear recoil band remains, for instance due to unaccounted for radioactive 
decays. If the signal in the nuclear recoil band modulates with the period and phase expected 
for a dark matter signal, however, it would be a strong indication of the nature of the interaction 
producing the signal. 

To date, only the DAMA (TJ [8] and CoGeNT [H [2] experiments have claimed an annual mod- 
ulation signal. The DAMA experiment, which uses target crystals of Nal(Tl), claims an 8.9a 
modulation with period 0.999 it 0.002 years and peaking at 146 it 7 days. The CoGeNT experi- 
ment, which uses a Ge target, has recently claimed an annual modulation signal with their first 
fifteen months of data. For energies ranging from 0.5-3.0 keVee,^ they observe a maximal modula- 



^ The notation keVee refers to the "electron equivalent energy in keV" , which is defined as the reconstructed recoil 
energy under the assumption that it is carried by an electron. For nuclear recoils, only part of the recoil energy 



tion with best-fit modulation fraction of 16.6 it 3.8%, period 347 it 29 days, and minimum at Oct. 
16ibl2 days. In this energy range, the significance is 2.8cr. When the CoGeNT data are fit with a 
dark matter signal in addition to a constant and exponential background, the best-fit dark matter 
mass is roughly 7-8 GeV with a ~ 10^^ pb, which is close to the region of parameter space that is 
consistent with DAMA [HEIE]. 

The light dark matter interpretation of CoGeNT and DAMA has been challenged by null results 
from other direct detection experiments, such as XENONIOO [lOj, XENONIO [ID IS], Simple [HI 
m, and CDMS [IS US]. The compatibility of the XENON and CoGeNT results has been discussed 
in greater detail in |17H21| . Reconciling CDMS and CoGeNT is more challenging because the 
two use the same target material and CDMS has reported an event rate significantly below that 
of CoGeNT in the same energy range; however, it has been claimed that errors in the energy 
calibration can potentially cause the discrepancy p2j. As we shall see, our conclusions on the 
modulation at CoGeNT will not depend strongly on these details. 

This work presents a detailed statistical analysis of the CoGeNT results using the publicly 
available data [23] . Section O introduces the statistical tests that will be used throughout the 



paper. Section III presents a model-independent analysis of the modulated and unmodulated rate. 



period, and phase. Section IV discusses the implications of these results for dark matter and shows 
that the hypothesis of a light, elastically-scattering WIMP (Weakly Interacting Massive Particle) 
is strained. In addition, the consistency of CoGeNT with CDMS-Ge, CDMS-Si, XENONIOO and 
DAMA is presented using an analysis that is independent of astrophysical uncertainties. 

II. COGENT DATA AND ANALYSIS TECHNIQUES 

We perform our analysis on the data that is available from the CoGeNT Collaboration upon 
request |23j . The first event was recorded on December 4, 2009 and the full data run spanned 
458 days, of which 442 were live. The known background in the energy region of interest arises 
from cosmogenic L-shell electron capture events. A complete description of the backgrounds, 
efficiencies, and experimental deadtime is included in Appendix [A} 

A binned chi-squared analysis of the data is done to confirm the 2.8cr significance for modulation 
in the 0.5-3.0 keVee energy bin [21 [3]. In addition, we use two alternate statistical tests that are 

is visible in the detector — an efTect that has to be corrected for by dividing the visible energy by a quenching 
factor — so that the energy threshold for nuclear recoils is higher than that for electron recoils. When referring to 
a nuclear recoil energy, we will use the notation "keVnr" . 



appropriate for searches of periodicities in data: the unbinned maximuni hkehhood method [23] and 
the Lomb-Scargle periodogram |25|l26j. Both these tests have been used, for example, in searches 
for periodicities in the solar neutrino flux from the Sudbury Neutrino Observatory [23 • Below, we 
review these statistical tests and note that, where the different techniques can be compared, they 
give qualitatively similar results. 

A. Binned Analysis 

We carry out a simple chi-squared analysis on events in a given energy range, binned in time. 
Each event, located at energy Ei, is reweighted by the efficiency at that energy, fcs{Ei)~^. The 
cosmogenic background contribution is determined for each bin and then subtracted from the 
number of events in the bin. A correction is applied to the time bins that overlap with the 
shutdown periods of the detector. The errors in each bin are treated as Gaussian,^ based on the 
original bin contents before reweighting or subtraction. 

The subtracted binned data are fit with a modulated spectrum of the form 

R{t) = Ao{l + Aicosiuj{t-tQ)), (1) 

where uj is the oscillation period, to is the phase, Ai is the modulation fraction, and ^o is the 
unmodulated rate. Note that the unmodulated rate may contain contributions from a dark matter 
component, as well as any other constant backgrounds. All times are taken relative to January 1^*, 
2010. We consider several types of fits, including: (1) all parameters are allowed to float, (2) the 
period is fixed to one year, ujq = 27r/year, (3) the period is fixed to ojq and the phase is fixed to 
to = 152 days, the value expected for dark matter in the standard halo model (SHM), and (4) the 
modulation fraction is set to zero (null hypothesis). 

B. Unbinned Analysis 

To maintain access to all the information in the time distribution of events, we carry out an 
unbinned maximum likelihood analysis. This method can be used to test the hypothesis that the 
excess data above backgrounds follow a rate distribution of the form in Eq. [T} If no assumption 



is made about the energy distribution of the data, as in Sect. Ill then it should be binned in 



^ This is a reasonable approximation for the number of time and energy bins used in this paper. 



energy, but not in time. The probability density function (PDF) for events in the energy range 
(Eiow; Ehigh = Elow + AE) is 

'■Ehieh 



m 



0.33 kg xAE/eff(E)i?(t)+ / '"''/cosmo(^,t)/cff(E) 



/gaps(t) , (2) 



where /eff(E) is the weighted average efficiency in the energy bin, fcosmoiE, t) is the model of the 
cosmogenic backgrounds, and /gaps(i) accounts for experimental deadtimes. See Appendix [A| for 



a more complete discussion of the background modeling. In Sect. IV, we carry out fully binned 
and unbinned analyses to establish whether the data is explained by dark matter. In this case, the 
PDF does not include the integral over energy. 
Using this PDF, the extended (log-) likelihood is 

f'tcnd 

2 log L{Ao, Ai,uj, to) = 2 ^ (l){ti) - 2 / " dt ^{t) , (3) 

j "^ ^start 

where tstart = —28 days and tend = 429.9 days for the data run in |2j, and the sum is over all data 
points in the sample. The quantity L{AQ,Ai,u},to) is maximized, subject to the constraints that 
Aq^i be positive and that < to < year. The significance for any particular hypothesis relative to 
any other is the difference between 2 log Lmax for each; this difference follows a x^ distribution. 

1. Lomb-Scargle Periodogram 

The weighted Lomb-Scargle technique is ideally suited to search for periodic signals in unevenly 
sampled data, such as that of CoGeNT. For data that are divided into A^ independent time bins 
with y(tj) data points each (z = 1, . . . , N), the Lomb-Scargle power for frequency / is given by 

„..._ 1 ( [EliW^{y{U)-y)cosu{U-T)]' [Eli W^y jU) - y) sm ujU - t)]' \ 
^^^ 2^2 1^ Y.li W. cos2 c.{U -r) ^ Ell W. sin^ u{U - r) J' ^ ' 

where y and a are the weighted mean and variance for the data in all the time bins and uj is the 
angular frequency. The phase factor r and weight factor Wi are given by 

, ,o ^ EliWism2ujti , ^^^ 1/af ._. 

tan(2^T) = -^.^— and W^ = -—^ , (5) 

2^i=iWi cos 2u}ti {^/(^i) 

respectively. Here, (Ji are the individual uncertainties in each bin. 

For a given energy range, the events are divided into eighty time bins of approximately six 

days each. In order for the Lomb-Scargle analysis to have a well-defined statistical interpretation, 

the contents of each bin must be large enough that the error on the number of events is well 

approximated by a Gaussian. As a result, each bin is required to contain ten or more events. 
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FIG. 1: Time-binned data in various energy ranges. Specifically (left) [0.5-1.5] keVee, (right) [1.5- 
3.1] keVcc. Overlaid are the best-fit to the modulation, as derived using the binned analysis, with 
free phase {solid red curve) and peak set at 152 days (dashed blue). The best-fit points correspond to 
Ao = 7.4 (7.5) events/day/kg/keVee, Ai = 0.14 (0.09) and to = 107 (152) days, for the phase free (to = 152 
days) for the lower bin and Aq = 2.7 (2.7) events/day/kg/keVee, Ai = 0.18 (0.14) and to = 116 (152) days 
for the higher. 

A simple algorithm is used to merge any bin that contains fewer than ten events with the next 
highest bin.^ This procedure is repeated until no bin has fewer than ten events. In addition, 
the centers of the bins are shifted to take into account any deadtime in the experiment. Finahy, 
the bin contents are efficiency-adjusted, the L-shell background in every bin is subtracted off, and 
the contents of the bin are converted to units of events/day /keVee. The error is based on the 
total (pre-subtraction) bin contents. This error is important for determining the weighting factors 
Wi. The power observed in the frequency uq = 27r/year can be converted to a significance for an 
oscihating signal. The probabihty of observing power P at any particular frequency in data that 
do not contain an oscillating signal is e~^, whereas the probability for observing power P at any 
frequency (including the appropriate trial factor) is approximately 1 — (1 — e~^)^, where A^ is the 
number of time bins j28j . 

III. A STUDY OF MODULATION 

The central goal of this work is to understand the properties of a potential modulation in the 
CoGeNT data. Therefore, we begin by applying the statistical techniques presented above to 
analyze the properties of the modulation, without any assumptions of its origin. We reproduce the 
results in [2], where a time-binned analysis is done in the energy ranges 0.5-0.9 keVee and 0.5- 



If the last bin has fewer than ten events, it is merged with the penultimate bin. 
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3.0 keVee. The results of [2] suggest that the region above 0.9 keVee exhibits a sizeable component 
of the modulation. For clarity, then, we divide the energy into two exclusive regions: a "low region" 
[0.5-1.5 keVee] and a "high region" [1.5-3.1 keVee], shown in Figurefl} 

The data clearly exhibit a modulation over a wide range of energies. The significance of adding a 
modulating term with a free phase (2 parameters), relative to the null hypothesis of no modulation, 
is Ax^ = 4.7 in the range [0.5-1.5] keVee, and 8.2 in the range [1.5-3.1] keVee. Furthermore, the 
improvement in Ax^ arising from adding a cosine with fixed phase (1 parameter) is 2.3 for the 
[0.5-1.5] keVee range, and 5.2 for [1.5-3.1] keVee. There is strong support for modulation in the 
high energy range, with little benefit in the low range. In addition, the high energy region prefers 
a phase that is different from that expected for a Maxwellian halo. 

The modulation in energies above 1.5 keVee is surprising, as the rate spectrum in this region 
had previously been interpreted as a constant background contribution [Ij. The unexpected nature 
of this modulation warrants a careful analysis of its properties and in the next two subsections, we 
apply additional tests to study its period, phase and amplitude. The last subsection presents the 
energy spectrum for the unmodulated and modulated rates as well as for the phase, assuming an 
oscillation period of a year. 
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FIG. 2: Significance of daily modulation in CoGeNT, as measured by the probability for the null (no 
modulation) hypothesis to give the observed amount of modulation, in a bin with energies ranging from 
Blow to Ehigh- Wc fit a model of the form (IT]) to the data after subtracting backgrounds and correcting 
for detection efficiencies and shutdown periods. The oscillation period w is kept fixed at (a) one solar day 
(24 hrs) and (b) one sidereal day (23.93 hrs); the average rate Ag, the modulation fraction Ai, and the 
phase to are free parameters in the fit. 



A. Oscillation Period 

The first step in characterizing the CoGeNT modulation is to determine the relevant time 
periods that show up in the data. The most obvious to check is evidence for daily modulation. 
While the daily modulation expected from dark matter is negligible in a detector like CoGeNT, 
many sources of background, such as those induced by radon decays and cosmic rays, can depend 
on the time of day. For instance, at night the atmosphere is colder and denser so that the secondary 
pions produced in cosmic ray interactions are more likely to lose energy through scattering before 
they decay [29] . 

Figure [2] shows the significance of modulation, under the assumption that the oscillation period 
is one solar day (24 hrs) and under the assumption that it is one sidereal day (23.93 hrs). The 
plots show results for different energy ranges i?iow through E'highi where E'low and -Ehigh range from 
0.5 and 3.0 keVee in steps of 0.1 keVee. Even though the best-fit solutions typically include about 
10-20% modulation, the statistical significance is very small, as indicated by the p-values in the 
plots. A few isolated energy intervals exhibit a modulation with more than 2a confidence, but 
the significance is much lower once the trial factor for this to happen anywhere in the considered 
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FIG. 3: Results of the Lomb-Scargle analysis. The left panel shows the Lomb-Scargle periodogram for the 
full energy range 0.5-3.1 keVee and the right panel plots the significance of annual modulation as a function 
of the considered energy range, from Ejow to Ehigh- 



energy range is included. We thus conclude that the CoGeNT data do not show evidence for 
diurnal modulation. 

More generally, we can also search for modulation with any frequency using the Lomb-Scargle 
technique. The results are shown in Fig. |3j The strongest modulation in the data has a period 
of one year; the Lomb-Scargle significance for annual modulation is above 3a if no trial factor is 
included, and around 90% with a trial factor. The right panel of Fig. |3] shows the significance of 
annual modulation, defined as the probability of obtaining the observed annual modulation from 
statistical fluctuations alone (not including a trial factor), as a function of the considered energy 
range from Ey^^ to i^high- There is no significant modulation below ~ 1.7 keVee, but the significance 
increases once higher energies are included. 

B. Phase and Amplitude for Annual Modulation 

Next, we consider the phase and amplitude of the modulation for a constant period of one year. 
An unbinned (in time) log-likelihood analysis is done for three energy ranges: low [0.5-1.5] keVee, 
high [1.5-3.1] keVee, and all [0.5-3.1] keVee. Figure [4] shows that the high energy data carry 
nearly the full weight of the analysis, and that the preferred phase is not Maxwellian, confirming 
the results of the binned analysis in Fig. [T] The modulation fraction in the high energy range is 
~ 20%, with a phase around 106 days. 

Figure [5] shows the significance of modulation over the null hypothesis in a range of energies 
from Flow to Ehigh, where Fiow, Ehigh each go from 0.5 to 3 keVee in steps of 0.1 keVee. The results 
for both the binned (left-hand column) and unbinned (right-hand column) analyses are shown to 
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FIG. 4: Likelihood analysis of the allowed regions in modulation and phase for different energy ranges: in 
[0.5-1.5] keVee {left), [1.5-3.1] keVee {middle), and [0.5-3.1] keVee {right). The contours are of Ax^ from 
the best-fit point, shown as •. 
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FIG. 5; Results of the modulation fit with lo — 27r/year, using both the binned and unbinned approaches. 
The upper plots allow the phase to float and the lower plots fix it to the value expected for the SHM 
(to = 152 days). The probabihty of the null (no modulation) hypothesis to fiuctuatc to the observed best-fit 
values is calculated from the Ax^ between the two best-fits, assuming 2 degrees of freedom for the upper 
plots and 1 for the lower. 



illustrate that the two methods are in very good agreement. The figure shows the significance when 
the phase is allowed to float in the upper row, and when the phase is fixed to Maxwellian in the 
lower row. 

The smallest p-values for the null hypothesis occur in the energy range 0.5-3.0 keVee. As in 
Fig. [3] (b), there is no significant modulation from Eiow = 0.5 to Ehigh ~ 1-7 keVee, but the signifi- 
cance starts to increase as Ehigh ^1-7 keVee. In the energy range where the modulation appears to 
be most significant, the phase and modulation fraction are both relatively stable; the best-fit phase 
falls consistently from 60-120 days, while the best-fit modulation fraction falls between 10-20%, 
for fits over the full energy range (0.5-3.0 keVee). 
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FIG. 6: Spectra for three different scenarios: (red square) the phase is allowed to float in the fit, (blue 
diamond) the phase is fixed to Maxwellian (152 days), and (green open circle) the phase is fixed to the 
best-fit phase (106 days) for a fit to the full data range 0.5-3.1 keVee. The spectra represent: total rate 
(left), modulation amplitude (middle) and phase (right). In all cases the error bars correspond to moving 
away from the best fit point by Ax^ — 1 (dark) or Ax^ = 4 (light), the horizontal offsets are for visualization 
purposes only. 

C. Spectra 



Having studied the basic properties of the modulation, we now consider the energy spectra of 
the unmodulated and modulated rate components, and the oscillation phase. The data are divided 
into four energy bins: [0.5-0.9], [0.9-1.5], [1.5-2.3] and [2.3-3.1] keVee. This separates the data 
into a bin dominantly below the cosmogenic peaks, one encompassing the cosmogenic background, 
and two evenly spaced bins above the cosmogenic peaks. For each bin, an unbinned (in time) log- 
likelihood approach is used to obtain the best-fit values for the unmodulated rate (^o); modulated 
amplitude {Aq * Ai), and phase (to)- 

Figure [6] shows the best-fit values as a function of energy, once the period is fixed to one year. 
Three different scenarios are considered: (blue) the phase is set to Maxwellian (to = 152 days), 
(green) the phase is fixed to the best-fit value for the fit over 0.5-3.1 keVee (to = 106 days), and 
(red) the phase is allowed to float bin-by-bin. Note that both the unmodulated rates and the 
modulation amplitudes do not differ dramatically between these three scenarios. The modulated 
amplitude has a large value (~ 1.5 cpd/kg/keVee) in the lowest energy bin, but very large error 
bars (±0.8 cpd/kg/keVee). The modulation amplitude flattens out considerably at higher energies 
around ~ 0.5 cpd/kg/keVee, with smaller error bars. When the phase is allowed to float, it remains 
relatively stable over the full energy range, with small variations around the best-fit value. The 
significance is highest in the last energy bin for the floating and best-fit phase scenarios; it is largest 
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FIG. 7: Same as Fig. [6J except for data from the high energy channel. 



in the second to last energy bin for the Maxwellian case. 

The fact that the modulation amplitude in Fig.[6]is non-zero even at large energies is a non-trivial 
and important feature of the spectrum. To see if this modulation amplitude eventually "turns off," 
we have also analyzed CoGeNT's high energy channel. The low and high energy channels in the 
data both measure pulse amplitudes in the range from 0.05 V and 0.25 V. However, the relation 
between the actual physical energy deposited and the measured voltage is different for the two 
channels. In particular, the response function at high energies is optimized to provide a good fit 
to the K-shell cosmogenics from ~ 4-12 keVee (see Appendix [A|) , but is not optimized for energies 
below ~3.2 keVee. Figure [7] shows the results of a log-likelihood analysis in equally-spaced bins 
above 4 keVee. The unmodulated rate spectrum indicates that there is an unexplained excess of 
events that is fairly constant with energy. In contrast to the low energy channel, the phase is highly 
unconstrained above ~4 keVee when it is allowed to float in the fitting procedure. This is due to 
the fact that there is no signifiant modulation in this energy regime, as illustrated by the spectrum 
of the modulation amplitude. 

Figure [6] shows that there is some modulation in the energy bin from 0.9-1.5 keVee, where the 
L-shell cosmogenic peaks are expected to dominate. As a check that the most dominant of these 
peaks, ^^Ge, is not oscillating, we plot the time-binned data for energies centered on this peak 
in Fig. M (red open circles). The black diamonds in the figure show the expected exponentially 
falling background, modeled using the procedure described in Appendix [^ The bottom panel in 
the figure shows the residual between the data and model, overlayed with the best-fit modulation 
in the 0.9-1.5 keVee energy bin from Fig. [6J As the energy range about E = 1.3 keVee is widened, 
the residuals come into better agreement with the blue line. In the narrow energy band, the 
statistics are too small to conclude whether there is a larger modulating component on the peak. 
As more data are collected and the statistics improve, it will be crucial to study the time variation 
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FIG. 8: Time variation for data centered on the ®^Ge L-shell peak for two different energy ranges. The top 
panel shows the predicted cosmogenic contribution using Eq. |A1| and the parameters given in Appendix [K\ 
(black diamonds), as well as the (efficiency corrected) time-binned distribution of the data (red open circles). 
A constant of 1.4 counts/day/keVee (see unmodulated spectrum in Fig.l6| has been added to the background. 
The bottom panel shows the residuals between the data and the model (red). The dashed blue line is the 
best-fit modulation in the range 0.9-1.5 keVee, obtained using the log-likelihood approach as in Fig. |6] 

of the ^^Ge line, especially as with more time it should decay away allowing easier access to any 
underlying modulation. 

IV. A DARK MATTER INTERPRETATION 
A. Dark Matter Fit for Standard Halo Parameters 



The results of the previous section suggest that the distribution of the modulated amplitude is 
not indicative of a conventional elastically-scattering WIMP. Furthermore, there is some discrep- 
ancy between the best-fit phase and that expected from the standard halo model (SHM). In this 
section, we will explore various dark matter (DM) models to see what best fits the observed data. 
Unless specified, the velocity distribution of the dark matter is Maxwell-Boltzmann with velocity 
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dispersion vq = 220 km/s and escape velocity Vesc = 550 km/s: 

f{v) oc (e-^'/^'o - e--'-/-o) e(t;esc - v) , (6) 

where u is the velocity in the galactic rest frame. More general velocity profiles will be considered 
in the following subsections. 

We carry out fits using an unbinned extended maxiniuni likelihood approach and a binned x^ 
analysis. For the unbinned method, we define a likelihood function that includes the dark matter 
signal, the cosmogenic backgrounds, and a constant background with floating normalization. The 
likelihood function accounts for efficiencies and shutdown periods. For the binned approach, the 
data is divided into five energy bins of equal size, spanning the range from 0.5-3.0 keVee. Within 
each energy bin, the events are partitioned in fifteen equal-sized time bins, each approximately one 
month wide. We subtract cosmogenic backgrounds and correct for the shutdown periods of the 
detector (the efficiencies are accounted for in the predicted dark matter signal). The error in each 
bin is based on the statistical uncertainty before background subtraction and deadtime correction, 
and is assumed to be approximately Gaussian, which is a reasonable assumption because no bin 
has fewer than six events. We carry out a x^ fit to these 75 bins, minimizing over dark matter 
parameters and systematic nuisance parameters. 

Several different signal and background scenarios are considered and the significance for each is 
summarized in Table |I} These results should be compared to a background-only fit where a constant 
rate is assumed in each of the five energy bins with no time variation. This fit gives x^ = 58.2 for 
70 degrees of freedom (d.o.f.). The separate scenarios are: 

• Spectrum -|- Modulation We attempt fitting elastic DM (eDM), varying a and m^, to 
both the energy spectrum and its time dependence. We assume an additional background 
contribution, constant in time and energy, and include its rate as a nuisance parameter, cq. 
The binned analysis gives x^ = 57.3 for 72 d.o.f.. The predicted and observed rates for 



the best-fit point obtained from the unbinned method is also shown in Fig. |9] and 10 The 
best-fit for the eDM scenario, which has a mass ~ 7 GeV, is marginally better than the 
time-independent background-only fit; the corresponding p- value is 0.64. 

Spectrum Only The second column of Table |l] shows an eDM scenario with a constant 
background, but this time fitting only to the unmodulated spectrum, i.e. ignoring time 
information. In this case, the binned least-squares analysis uses 50 bins and gives x^ = 50.8 
for 47 d.o.f. 
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Scenario Spect+Mod Spect only Mod only Mod only iDM 

(c, > 0) Mod only 

d.o.f. 72 [n/a] 47 [n/a] 

cr/lO-*! cm^ 13.8 [8.9] 10.1 [8.2] 
m^/GeV 7.2 [8.1] 7.7 [8.2] 

(5/keV 



c,/(cpd/kg/keVee) 2.5 [2.5] 2.5 [2.6] 



68 




68 




67 


6.0 8.6 64 


10.0 12.0 16.3 


24 


('■'] 




(''■'] 




('■'] 







-4.6 




0.4 


1.5 




-1.0 




0.4 


2.3 




1.1 




1.2 


2.3 




1.7 




1.5 



57.3 [n/a] 50.8 [n/a] 53.7 51.4 



51.3 



TABLE I: The best-fit dark matter parameters for the binned [unbinned] analyses of various DM scenarios. 
The numbers should be compared to a background-only fit that gives x^ — 58.2 for the 70 degrees of freedom 
(see text for details). 

• Modulation Only The third and fourth columns of Table|I]only fit to the modulation in the 
data. The fits are done for an eDM plus background hypothesis, where the time-independent 
background is allowed to vary freely in each of the five energy bins (i.e. , there are five nuisance 
parameters Cj). For the case where these constant contributions, presumably coming from 
some unidentified background, are restricted to be physical (cj > 0), the fit gives x^ = 53.7 
for 68 d.o.f. If instead they are allowed to float, x^ = 51.4, again with 68 d.o.f. Notice that 
when fitting the modulation, the best-fit dark matter mass is ~ 10 — 12 GeV. 

Table |I] shows that although the data are consistent with the SHM DM hypothesis, the inclusion 
of such a DM component does not greatly improve the fit compared to the hypothesis of a back- 
ground that is constant in time. Furthermore, there is slight tension in the DM interpretation — the 
modulation in the data by itself favors heavier dark matter masses than the spectrum. This be- 



havior is also illustrated in Figure 11 where we show the preferred region in DM mass m-y and 



elastic spin-independent scattering cross section crsi for various fits to the CoGeNT data. The plot 
confirms the slight (but not statistically significant) tension between the DM masses preferred by 
the energy spectrum observed in CoGeNT and the annual modulation. A fit to the modulation 
data alone can exclude the hypothesis of no DM at low confidence level (light red contours), how- 
ever this requires an unphysical background (i.e., Cj < 0). If positive background is required, the 
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FIG. 9: Comparison of the CoGeNT data to the predicted event spectrum for elastically (sohd dark red) 
and inelasticahy (light red dashed) scattering dark matter. In both cases, results for the best-fit dark matter 
parameter are shown. In the elastic case, the fit is done using the unbinned maximum likelihood approach, 
including energy and timing information for each event as well as a constant background. For the inelastic 
case, a binned analysis is done with a background that can float in each of five energy bins. (For the iDM 
case, the signal is shown, but not the fitted constant background.) 



modulation provides only an upper bound. We also compare the CoGeNT-preferred regions to 
the exclusion limits from XENONIOO [TO] and GDMS [16], and confirm the well-known tension 
between these results. To compute the DM parameter region favored by DAMA [7], we follow the 
procedure described in [30j, assuming no channeling |31tl32| and assigning a 10% systematic uncer- 
tainty to the DAMA quenching factors [331^35] . This systematic uncertainty has a strong effect on 
the horizontal extent of the DAMA region. For standard halo parameters, the DM interpretations 
of DAMA and CoGeNT are inconsistent. To some extent, this tension can be relaxed for halo 



parameters different from the ones we chose in Fig. 11 |361 I37j . In addition, one might speculate 
that systematic uncertainties in either DAMA or CoGeNT can be larger than what was assumed 
here. For instance, if the low-energy event excess in CoGeNT is only partly due to dark matter 
and partly due to some other source, higher DM masses may become allowed. 

Because an elastically scattering WIMP with a Maxwell-Boltzmann velocity distribution only 
produces a large modulation amplitude at low energies (< 1.5 keVee), it is worthwhile to consider 
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FIG. 10: Comparison of the CoGcNT data to the predicted time-dependent event rate for elasticaUy (sohd 
dark red) and inelasticahy (hght red dashed) scattering dark matter, using the same fitting procedures as 
for Fig. [9| 



non-standard WIMP scenarios as an explanation of the CoGeNT modulation. Inelastic dark matter 
(iDM) is an example of a DM scenario with increased modulation and a preference for events at 
high recoil energy |38] . We can repeat the exercise above under the assumption that DM scatters 
inelastically |39j . However, there is no preference for iDM from fits to the unmodulated spectrum, 
which has many events below where an inelastic contribution is expected - see Fig. |9j There is, 
however, a preference for iDM if one fits to only the modulation (see last column in Table ^. 
Interestingly, for the iDM hypothesis, there is no improvement in the fit if the background is 
allowed to be unphysical (cj < 0). Instead, the iDM model enables a similar modulation spectrum 
as the eDM fit, but with physical backgrounds. For the modulation-only fit, iDM gives x^ = 51-3 
for 67 d.o.f. The spectra for this best-fit point are shown in Fig. [TUJ clearly, iDM can produce 
modulation in all energy bins, but at the cost of explaining only part of the event excess observed 
by CoGeNT. We note that inelastic WIMPs are highly sensitive to non-Maxwellian properties. If 
the phase is truly shifted away from the Maxwellian 152 day peak, this could be indicative of such 
halo properties. The CoGeNT best-fit phase is mid-April, during the time when XENONIOO had 
increased levels of noise; therefore, a signal localized in this time could have been missed. Because 
iDM in the presence of a stream can lead to narrow peaks for brief periods of the year |40j , further 
running of XENONIOO in April should clarify this situation. 
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FIG. 11: Preferred regions and exclusion limits at 90% and Scr confidence level in the m^-a plane for spin- 
independent dark matter-nucleon scattering assuming a standard Maxwell-Boltzmann halo with escape 
velocity Uosc = 550 km/s and velocity dispersion va = 220 km/s. Filled red (dark gray in B/W) contours are 
obtained from an unbinned maximum likelihood fit to the CoGeNT data, using both the energy and timing 
information for each event. (A fit using only energy information gives practically identical results.) The 
unfilled red (gray) contours are from a binned x^ analysis, using only the timing information and leaving the 
energy spectrum completely unconstrained (light red/light gray contours) , or requiring the predicted energy 
spectrum to remain below the observed one (dark red/dark gray exlcusion limits). The orange (light gray) 
region shows the masses and cross sections preferred by DAMA [7] if the quenching factors are assigned 
a 10% uncertainty [5D1 |33J [3S], and the blue and green contours indicate the 90% exclusion limits from 
CDMS [16] and XENONIOO ^U\, respectively. 

B. Varying the Halo Parameters 



Next, we explore whether the CoGeNT data are compatible with a general class of equilib- 
rium velocity distributions that extend beyond Maxwell-Boltzmann. In particular, we consider 
distributions of the form 



f{v) oc (e-^'/^o _ e--?BcK)fc e(„^^^ _ ^) , 



(7) 



where A; is a power-law index, Ugsc is the escape velocity, and vq is the dispersion. Note that k = 1 
is just the Maxwell-Boltzmann- like halo. This velocity distribution models the behavior of double 
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FIG. 12: Modulation amplitude in the range 1.5-3.1 keVee as a function of dark matter mass, where the 
dark matter cross section is normalized to fit the modulation amplitude in the first bin (left) and over the 
whole energy range (right). The colors indicate different spectral indices for Eq. It] k=l (blue), ^—2 (pink), 
k=3 (green). The regions between (above) the solid (dashed) lines indicate points that overpredict the 
unmodulated rate by at least 2a from 0.5-1.5 keVee (1.5-3.1 keVee). The solid colored bands are the only 
regions consistent with the unmodulated rate spectrum. The gray band is the modulated amplitude with 
Ict error bars for the 1.5-3.1 keVee region. 



power-law density profiles and corresponds to results found in high-resolution simulations of the 
Galactic halo, when A; ~ 2 [37j. The fact that simulations support a power-law index greater 
than one suggests that the number of high velocity particles on the tail of the distribution may be 
suppressed relative to the expectation for Maxwell-Boltzmann halos. 

To study how the CoGeNT predictions are affected by variations in the halo parameters, distri- 
butions with k = 1, 2, 3 are considered and a random scan is done with Wcsc £ [500, 600] km/s |41j 
and vq G [180, 280] km/s |42H44j . For each randomly selected set of halo parameters and dark mat- 
ter mass, the modulated and unmodulated rates are evaluated in two energy bins: 0.5-1.5 keVee 



and 1.5 - 3.1 keVee. The results are summarized in Fig. 12 for the case where the cross section is 



normalized to give the measured modulation amplitude in the first bin (left) and the case where it 
is normalized to the total amplitude from 0.5-3.1 keVee (right). The three colors represent three 
different spectral indices, with k=l Maxwellian-type in blue, k=2 in pink, and k=3 in green. The 
gray band is the best-fit modulation amplitude (±1(7 error bars) for 1.5-3.1 keVee, obtained using 
the log-likelihood method. 

The data support a modulation amplitude of ~ 0.38 it 0.16 cpd/kg/keVee in the high energy 
bin. For the elastic scattering case considered here, only dark matter masses greater than ~ 9 
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GeV for k=l and ~ 10.5 GeV for k=3 yield a modulation amplitude within a standard deviation 
of the measured value. While a heavier dark matter mass increases the modulation amplitude at 
high energies, it also increases the unmodulated rate, making it conflict with the rates measured 
by CoGeNT. For the wide range of halo and dark matter parameters considered here, only masses 
less than ~ 7 — 9 GeV are consistent with CoGeNT's unmodulated spectrum. Unfortunately, none 
of these points give a sufficient modulation at high energies. The results shown here emphasize 
the underlying tension in a dark matter interpretation of the CoGeNT data: namely, one must 
explain both an excess in the unmodulated rate below ~ 0.9 keVee and a significant modulation 



above ~ 1.7 keVee. As Fig. 12 highlights, a dark matter candidate scattering elastically off an 
equilibrated isotropic halo cannot satisfy both requirements. However, this does not mean that it 
is impossible. In particular, a two-component halo that is described by a Maxwell-Boltzmann at 
low velocities and a stream at high velocities might allow consistency. A modulation phase that is 
significantly different from 152 days would further point to a more exotic halo model. 

C. Model-Independent Comparisons 

In this subsection, we explore the constraints from other experiments on the CoGeNT mod- 
ulation, assuming it arises from elastic dark matter. Comparing rates between different direct 
detection experiments with different target nuclei is non-trivial because each probes a different 
range of dark matter velocities. However, a means of comparing the results of different experi- 
ments independent of halo models has recently been proposed [IS].^ For elastic spin-independent 
scattering, a signal in the range [E[^^, EJ^-^^j^] at Experiment 1 arises in Experiment 2 in the energy 
range 

rp(2) -p(2) , _ /ip4^r-p{l) -p(l) 1 ff.. 

L-^low'-^highJ ~ 2 »^(2) L low' -^highJ ' \°) 

Li^ IVlrp 

(i) 

where M^ is the mass of the target nucleus in each experiment and fn is the DM-nucleus reduced 
mass for each experiment. For a rate, dRi/dEn, observed at Experiment 1, the rate expected at 
Experiment 2 is 



dEn ' '' 41) 2 (^\<L^ \ dER \^2 ^a) - 

1 VmK' V 



For related work see 1461, 
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Bin 


CoGeNT 


Ge 


Na (Q=0.3) 


Si 





Xc 


1 


[0.5,0.9] 
0.90 ±0.72 


[2.3,3.8] 
0.23 ±0.18 


[1.5,2.5] 
0.078 ±0.062 


[4.5,7.6] 
0.035 ±0.028 


[5.8,9.9] 
0.011 ±0.009 


[1.4,2.3] 
0.72 ±0.58 


2 


[0.9,1.5] 
0.37 ±0.55 


[3.8,6.1] 
0.1 ±0.149 


[2.5,4.0] 
0.035 ±0.052 


[7.6,11.9] 
0.015 ±0.023 


[9.9,15.6] 
0.005 ± 0.008 


[2.3,3.7] 
0.31 ±0.46 


3 


[1.5,2.3] 
0.48 ±0.22 


[6.1,8.9] 
0.136 ±0.063 


[4.0,5.8] 
0.049 ±0.022 


[11.9,17.5] 
0.021 ±0.01 


[15.6,22.8] 
0.007 ±0.003 


[3.7,5.4] 
0.41 ±0.19 


4 


[2.3,3.1] 
0.27 ±0.23 


[8.9,11.6] 
0.08 ±0.068 


[5.8,7.6] 
0.029 ±0.025 


[17.5,22.8] 
0.013 ±0.011 


[22.8,29.8] 
0.004 ± 0.004 


[5.4,7] 
0.23 ±0.2 



TABLE II: Predicted modulation amplitudes for example nuclear targets, given the best-fit values for Co- 
GeNT assuming a Maxwellian phase. The units are in counts/day/kg/keVnr for all columns, except those 
labelled CoGeNT and Na where they are counts/day /kg/keVee; in the case of Na a quenching factor of 
Q=0.3 has been applied. The equivalent energy ranges and rates for other targets are shown, assuming 
m^ — 7 GeV and spin-independent scattering cross sections proportional to A^. Note that we have not 
included detector efficiencies or mass fractions in any of the predicted rates. 



Here, 



a 



(0 



K»f/pZ« + /„(A«-Z«) 



(10) 



where k is the mass fraction for the target element in question and fpfn) is the coupling strength 
of dark matter to protons (neutrons). Fi is the nuclear form factor for each experiment. 

Tables |II] and |III| show the ranges of energies at other experiments that correspond to the Co- 
GeNT energy bins: [0.5, 0.9], [0.9, 1.5], [1.5, 2.3], and [2.3, 3.1] keVee. Note that these energies 
are given in "electron equivalent" and correspond to [2.3, 3.8], [3.8, 6.1], [6.1, 8.9], and [8.9, 11.6] 
in nuclear recoil energies. These tables also show how the CoGeNT modulation amplitude in each 
energy bin translates to other experiments, assuming a 7 GeV WIMP with spin-independent scat- 
tering proportional to A^. (Note that we have not included detector efficiencies or mass fractions 
in any of the predicted rates.) Let us consider each experiment in turn. 

CDMS-Ge: A direct comparison can be made between the CoGeNT and CDMS count rates 
because they both have germanium targets. Using the results of the low-energy analysis of the 
CDMS experiment |16j . we calculate an upper limit for the rate in each detector such that it has a 
1.3% probability of having a lower rate. This gives a probability of 10% that any one of CDMS's 
eight detectors has a lower rate than is observed. In each of the five energy bins, the strongest 
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Bin 


CoGeNT 


Ge 


Na (Q=0.3) 


Si 





Xe 


1 


[0.5,0.9] 
1.4 ±0.79 


[2.3,3.8] 
0.36 ±0.2 


[1.5,2.5] 
0.12 ±0.07 


[4.5,7.6] 
0.054 ±0.03 


[5.8,9.9] 
0.018 ±0.01 


[1.4,2.3] 
1.1 ±0.6 


2 


[0.9,1.5] 
0.84 ±0.59 


[3.8,6.1] 
0.23 ±0.16 


[2.5,4.0] 
0.079 ±0.055 


[7.6,11.9] 
0.035 ±0.024 


[9.9,15.6] 
0.012 ±0.008 


[2.3,3.7] 
0.70 ±0.49 


3 


[1.5,2.3] 
0.46 ±0.24 


[6.1,8.9] 
0.13 ±0.068 


[4.0,5.8] 
0.047 ±0.024 


[11.9,17.5] 
0.021 ±0.011 


[15.6,22.8] 
0.007 ±0.004 


[3.7,5.4] 
0.39 ±0.21 


4 


[2.3,3.1] 
0.66 ±0.24 


[8.9,11.6] 
0.20 ±0.07 


[5.8,7.6] 
0.072 ±0.026 


[17.5,22.8] 
0.032 ±0.011 


[22.8,29.8] 
0.011 ±0.004 


[5.4,7] 
0.57 ±0.21 



TABLE III: Same as Table ITTl except assuming a best-fit overall phase of 106 days. 



limit from all the detectors is chosen and we treat this as a 90% confidence limit. ^ Figure 13 
shows that the count rates at CDMS are not low enough to constrain the CoGeNT modulation. 
However, the count rates are low enough that there should be modulation at a very high level in 
CDMS. Thus, even weak modulation constraints from CDMS could be very powerful; conversely, 
modulation should be apparent in a dedicated analysis, even with existing count rates. 

CDMS-Si: CoGeNT's modulation above 1.5 keVee should appear above threshold at CDMS- 
Si. Using the results from [47j, with 88 kg-days and assuming an efficiency of 0.2, CoGeNT's 
observed modulation rates map into a minimum of 3.3 ± 1.4 (5.0 ± 1.5) events for a Maxwellian 
(106 day) phase. CDMS-Si sees no events below 50 keVnr, so this predicted rate is borderline, 
but not excluded. Two further conclusions can be made. First, any signal should essentially 
be 100% modulated (i.e., no sizeable constant piece). Second, invoking significant interference 
between proton and neutron couplings [481 H9] to reconcile CoGeNT with XENON is untenable. 
Taking fn = — 0.7 fp turns the relative boost between Si and Ge (the ratio of Cy's) all the way up 
to 11, putting any observable modulation in the 1.5-3.1 range in serious conflict with the results 
form CDMS-Si. Note, that the low energy calibration of the silicon detectors may be subject to 
corrections [M], which would shift the energy threshold at CDMS-Si to higher recoil energies. We 
have taken a fairly conservative approach in comparing CoGeNT to CDMS-Si, for instance only 
considering events above 11.9 keV in CDMS, and taking the efficiency to be a flat 20%, and we do 
not expect the above tension to be greatly alleviated by potential threshold corrections. 

CRESST-O: CoGeNT's signal in the [0.5-1.5] keVee region translates to an energy range that 



^ The probability that the particular detector that sets the limit has a strong downward fluctuation is small, and so 
the confidence is actually better than 90%, but we treat it as a 90% C.L. to be conservative. 
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FIG. 13: Upper limits from CDMS (red) compared with modulation rates from CoGeNT assuming a 
Maxwellian phase {blue diamond) and the overall best-fit phase (green circle). 



falls primarily below the threshold of many of CRESST's individual detectors. However, the [1.5- 
3.1] keVee range should easily appear at CRESST above 15 keVnr in the oxygen band. With 600 
kg-days of exposure and an efficiency of 80%, a modulated signal of 8.75 it 3.77 (13.1 it 4.0) events 
is expected at CRESST, if the CoGeNT signal is due to a light elastic DM. This is consistent with 
the ~ 30 events that have been observed at CRESST [50j. No significant modulation has yet been 
reported. 

XENONIOO: The energy range for CoGeNT's high-energy modulation is relevant for 
XENONIOO. In particular, 1.5 keVee corresponds to 3.7 keVnr in XENONIOO (again, for a 7 
GeV WIMP), where the scintillation efficiency has been studied. Taking an exposure of 1450 kg- 
days (48 kg X 100.9 days x overall efficiency of 0.3), we predict 1020 it 470 (970 ±510) events from 
the third bin and 550 it 470 (1350 it 490) events from the fourth bin before taking into account the 
SI cut. 

Both of these bins are below the SI threshold for XENONIOO of four photo-electrons (PE), so 
we must calculate the probability that an upward fluctuation would occur. Both these bins occur 
above 3 keVnr, for which results are available for Lcff. Taking a value of Lefr=0.07 (approximately 
the lower boundary as measured by |51j). we predict efficiencies of 0.015 and 0.05 for the third and 
fourth bins, respectively. With these efficiencies, one would have expected 15.4 it 7 (14.7 it 7.7) 
and 27.4 it 23.3 (66.9 it 24.3) events based on CoGeNT's modulation in the third and fourth bins, 
respectively. Reducing the Poisson efficiencies to an acceptable level requires that Lgff^ 0.05 and 
^ 0.04 for the third and fourth bin, which lower the rates by a factor of approximately 3 and 
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6, respectively. Whether such small values of Leg are possible or not is still a subject of active 
discussion in the literature [T71 - [2n [5TH57] 

DAMA: Finally, we compare the measured modulation spectrum at DAMA with that of Co- 
GeNT. For a quenching factor of QNa = 0.3 and m^ ~ 7GeV, the CoGeNT modulation energy 
range corresponds roughly with DAMA's (see Tables II|III). The total modulation observed in the 



CoGeNT range 0.5-3.1 keVee yields a modulated rate of 0.04 ± 0.017 cpd/kg (0.065 ± 0.018) for a 
Maxwellian (best-fit-106 day) phase, which compares with 0.0444 it 0.0052 cpd/kg (for both MW 
and best-fit 146 day phases) at DAMA, assuming a conventional spin-independent mapping. While 
the total modulated rate at CoGeNT is roughly consistent with DAMA's, the energy ranges at 
the two experiments do not line up exactly and the predicted rate at DAMA based on CoGeNT's 
signal is somewhat smaller than what is observed independent of the astrophysical model. This 



is illustrated in Fig. 14 since the energy bins are comparable in size to the energy smearing at 
DAMA we ignore the effect of smearing. As previously noted |481 H9] , taking the proton and neu- 
tron spin-independent couplings to interfere can favor light targets and correct this. Such effects 
can happen through interactions via heavy fermionic mediators |49j or through Z's [58]. However, 
taking /„ = — 0.7/p to maximize the suppression at XENON boosts the modulation at DAMA 
relative to CoGeNT by a factor of six relative to the case of fn = fp, which seems in conflict with 
the data in hand. 

We also consider the specific scenario of [31], with a 7 GeV WIMP and a large quench factor 
in sodium, Q^a = 0.5 [34J. In this case, the lowest two bins at CoGeNT map roughly into the 
DAMA range 2-6 keVee, yielding 0.019 ± 0.015 (0.036 ± 0.016) cpd/kg for a Maxwellian (106 day) 
phase. The higher energy bins (1.5-3.1 keVee) map into the range of 6.6-12.7 keVee at DAMA, 
and predict a modulation of 0.02 it 0.0085 (0.03 it 0.009) cpd/kg to be compared with the observed 



—0.0008 lb 0.0064 cpd/kg observed in the 6-14 keVee range, see Fig. 14 Consequently, this mass 
and quenching factor are in tension with the high energy data. While this confiict is at the 2 (2.7) 
a level for Maxwellian (106 day) phase, ignoring it requires one to essentially ignore modulation 
that is as significant as the modulation that one is taking seriously. As a result, we believe that 
this particular scenario does not give a good fit to the data. 

1. Summary of Halo- Independent Conriparisons 

To summarize, this subsection compares results between different experiments in a manner that 
is independent of astrophysical uncertainties. The most direct comparison is between CoGeNT 
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FIG. 14: Comparison of CoGeNT and DAMA modulation amplitudes. The comparison is independent of 
astrophysics but assumes 7 GeV DM that scatters elastically and a quench factor in sodium of 0.3 (left plot) 
and 0.5 (right plot), since the energy bins are comparable to the DAMA resolution, we ignore the effects 
of energy smearing at DAMA. In both plots, the blue diamonds (open green circles) denote the CoGeNT 
prediction for the modulation spectrum at DAMA, assuming a Maxwellian (best-fit) phase. The black (gray) 
points are the results of a two bin (multi-bin) analysis by DAMA [5] . 



and CDMS-Ge; the results of the two experiments are consistent only if CDMS's rate is modulated 
at nearly 100%. Such a modulation should be easily visible in the CDMS data. 

Ultimately, while there is rough agreement between the size of the CoGeNT modulation and 
the DAMA modulation, the energy range over which the modulation is spread is in conflict with 
previous interpretations jMj invoking a large sodium quenching factor, because this disregards the 
modulation at high energies, which is as statistically significant as that in the lower energy range. 

The presence of modulation in the high energy range results in the greatest tensions with other 
experiments. The absence of a signal at CDMS-Si requires the signal to be highly modulated, while 
XENONIOO should have seen a signal unless Lefj is significantly smaller than the measurements 
of [51]. Invoking interference between protons and neutrons to alleviate XENONIOO constraints 
exacerbates tensions with CDMS-Si. 

The comparisons made in this subsection are only valid in the context of spin-independent 
scattering where the rate is proportional to A^. We have seen that any interpretation in this 
context is challenging because of the high-energy modulation in CoGeNT. Other models, such as 
spin-dependent or inelastic interactions, fall outside the scope of this analysis. 
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V. CONCLUSIONS 

We have studied the event rate from the CoGeNT experiment over the first 458 days of running. 
CoGeNT continues to see an unexplained excess at low energies and now claims evidence for a 
modulating spectrum. The presence of modulation makes the possibility of a WIMP explanation 
more compelling, but also requires a serious and thorough discussion of its implications for theory 
and other direct detection experiments. 

In this work, we confirm the original findings of [2j of a modulation in the low-energy (0.5-3.1 
keVee) data of CoGeNT, with a significance of 99% or 2.6cr (99.7% or 3a) for a Maxwell-Boltzmann 
(best-fit) phase. However, the details of the modulation make it somewhat confusing. The signifi- 
cance for modulation in the low energy range (0.5-1.5 keVee), where the signal from a light WIMP 
would be present, is only 87% (90%) for a Maxwell-Boltzmann (best-fit) phase. In contrast, the 
significance is 97.7% (98.3%) in the higher energy range (1.5-3.1 keVee). The absence of significant 
modulation below ~ 1.5 keVee is not troubling on its own, because modulation fractions can be 
larger at higher energies, and cosmogenic backgrounds at 0.9-1.5 keVee are significant. Nonethe- 
less, one cannot claim that the modulation in the low-energy regime gives evidence for a model 
while simultaneously disregarding the much more significant modulation in the high-energy regime 
— i.e., if a model does not explain a significant part of the high energy modulation, it cannot be 
claimed to explain the modulation. The significance of modulation in the 1.5-3.1 keVee region 
is confirmed by a variety of statistical techniques, and is an important contributor to the total 
modulation significance in the 0.5-3.1 keVee energy range. 

Many elastic dark matter models with Maxwellian halos do not give any significant modulation 
above 1.5 keVee for masses below ~ 7 — 9 GeV. We have performed a broad study of elastic 
WIMP scenarios with Maxwellian and NFW-consistent halos and have found no models that give 
modulation at high energies, while not exceeding the unmodulated count rate at lower energies. 
Halo models with debris flows [^ or streams might allow for this, however. Our best-fit WIMP 
point with a Maxwell-Boltzmann halo does not significantly improve the fit compared to a constant 
background. As a result, we believe that attempts to understand the CoGeNT modulation with a 
Maxwellian halo are likely to be unreliable and lead to erroneous conclusions. 

To that end, we have attempted to employ techniques that are independent of the halo model 
|45) when comparing the CoGeNT results to other experiments. We find that a direct comparison 
to CDMS-Ge allows the modulation, but predicts that a significant modulation should appear at 
CDMS as well. Note that small errors in energy, while important for interpreting a rapidly falling 
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background, should not affect a general modulation analysis such as this one. 

Other experiments also provide insight on the modulation in the 1.5-3.1 keVee range. In par- 
ticular, the absence of a signal at CDMS-Si suggests that any signal in this range should be ~ 
100% modulated. A study of XENONIOO is intriguing, because the modulation appears in a range 
that has been more directly calibrated with measurements of Lefj. We find that any significant 
modulation in the 1.5-3.1 range should have shown up at the 10-30+ event level at XENONIOO, 
unless the value of Lgfj is significantly lower than what has been found by [51] . 

A halo model-independent comparison to DAMA shows generally good agreement between the 
event rates for a light WIMP and Q^a ~ 0.3. However, the modulation at CoGeNT is in conflict 
with DAMA if one assumes QNa = 0.5, disfavoring previous interpretations utilizing a Maxwellian 
halo. 

Inelastic dark matter models can provide a good fit to the modulation, but not the exponentially 
falling rate at low energies. iDM in the presence of stream(s) or debris flows could lead to the 
highly-modulated signal that is observed. If the signal is only present in April, it could explain 
why XENONIOO did not observe a significant rate, as data taking in that period was limited by 
excess noise. 

In summary, the modulation at CoGeNT is an intriguing piece of the puzzle, but raises as 
many questions through the high-energy modulation as it answers. Future data from CoGeNT is 
imperative for understanding the presence or absence of modulation at low and high energies and 
thus, shedding light on whether~and what type of-dark matter could give rise to it. 
Note Added At the time of the completion of the manuscript, the authors became aware of related 
papers j60| I61j . which explore implications of the CoGeNT data for particle and astrophysical 
models. The authors of |60| carried out an analysis on the two bin data presented in [2j, while |61) 
use the full data set analysed in this work. Our conclusions are in qualitative agreement. 

After this manuscript had been submitted to the e-print arXiv, the CoGeNT collaboration 
released new preliminary results which indicate that part of the low-energy excess events may be 
due to surface backgrounds [62]. If confirmed, this would imply that the parameter region favored 
by the CoGeNT event spectrum would become larger and move to lower cross sections 
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Appendix A: CoGeNT data and background 



The known background in the energy region of interest arises from cosmogenic electron capture 
events, which can be modeled as a sum of decaying gaussians. A radioactive element, with Gaussian 
peak at £"0, width a, and half- life t^^ has an energy spectrum of 



/peaks (,-C!/ ; '•j — -''atoms 



log22_t/ti/2 1 ^-[E-Eo)/2a^ 



(Al) 



ti/2 v2vr(T 

Here, A^'atoms is the total number of atoms expected to decay in the detector. The number of 
K-shell decays is obtained by fitting to the peaks in the high energy spectrum. These numbers are 



provided in the public release of the CoGeNT data and are summarized in Table IV, where they 
have been corrected for efficiencies. We have been able to reproduce A'^atoms for the ®^Ge peak at 







L-shcU 


K-shell 


Isotope 


ti/2 (days) 


-^^atoms 


Ef) (keVcc) 


a (keVee) 


-^^atoms 


Ef) (kcVcc) 


a (kcVec) 


^3 As 


80 


14.7 


1.41 


0.0777 


133 


11.1 


0.120 


68Ge 


271 


736 


1.30 


0.0770 


6460 


10.4 


0.117 


esQa 


271 


60.9 


1.19 


0.0764 


553 


9.66 


0.114 


«^5Zn 


244 


243 


1.10 


0.0759 


2250 


8.98 


0.112 


56Ni 


5.9 


1.78 


0.926 


0.0749 


17.2 


7.71 


0.107 


56,58Co 


71 


10.9 


0.846 


0.0744 


107 


7.11 


0.104 


"Co 


271 


2.98 


0.846 


0.0744 


29.3 


7.11 


0.104 


55pg 


996 


51.8 


0.769 


0.0740 


489 


6.54 


0.102 


54Mn 


312 


24.3 


0.695 


0.0736 


238 


5.99 


0.100 


5icr 


28 


3.38 


0.628 


0.0732 


33.5 


5.46 


0.0975 


49 p. 


330 


17.2 


0.564 


0.0728 


172 


4.97 


0.0953 



TABLE IV: Data used to model the cosmogenic background for both L-shcll and K-shell decays. Vatoms is 
the number of atoms in the detector expected to decay via each mode (before efficiency corrections) , E^ is 
the binding energy, a is the energy resolution, and ti/2 is the half-life for the decay. 
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10.4 keVee using our own fitting procedure. The binding energies for tlie K-shell peaks are given 
in [1] and tlie resolution is given by tlie formula in |65j . with parameters from [1]. 

The expected number of L-shell decays at lower energies is related to the number of K-shell 



decays |66j : therefore, the third column in Table IV can be obtained from the measured quantities 



in the sixth column. The L-shell cosmogenics are relevant from ~ 0.5-1.7 keVee, where the dark 
matter signal is expected to dominate. The dominant contribution in the first year of running 
comes from ^^Ge, followed by ®^Zn. 

The observed energy spectrum at CoGeNT is affected by an efficiency that captures the decreas- 
ing sensitivity to signal near the experiment's energy threshold. The efficiency is approximately 
flat at ~ 0.87 between 0.7 - 3.0 keVee, and then drops to 0.75 by 0.5 keVee. Above 4 keVee, it is 
0.94. In this work, a spline interpolation of the efficiency data points is used, denoted by /efr(-E'). 

The CoGeNT data run spanned 458 days, of which 442 were live. The time gaps in which no 
data were taken must be properly accounted for in a study of annual modulation. To account for 
these outages, the following function is introduced: 

/gaps(t) = (0(67 -t) + e{t - 74)) (9(101 -t) + e{t - 107)) (9(305 - t) + 9(i - 308)) , (A2) 

where t is measured in days. Therefore, the complete expression for the spectrum of cosmogenic 
peaks is 

fcosmo{E,t) = fpea.ks{E,t)fgaps{t)fcs{E). (A3) 
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